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ABSTRACT 

We report on the detection of the two shortest period non-interacting 
white dwarf binary systems. These systems, SDSS J143633. 29+501026. 8 and 
SDSS J105353. 89+520031.0, were identified by searching for radial velocity variations 
in the individual exposures that make up the published spectra from the Sloan Digital 
Sky Survey. We followed up these systems with time series spectroscopy to measure 
the period and mass ratios of these systems. Although we only place a lower bound 
on the companion masses, we argue that they must also be white dwarf stars. With 
periods of approximately 1 hour, we estimate that the systems will merge in less than 
100 Myr, but the merger product will likely not be massive enough to result in a Type 
la supernova. 

Subject headings: white dwarfs — binaries: close, spectroscopic 



1. Introduction 



are the end 
o of all stars 



White dwarf stars (WDs 
point of stellar ev olution for 
(jWeidemann1l2000l ) and store the archaeologi- 
cal record of the Galaxy. WDs in binaries are 
particularly rich systems to study. Because of 
their intrinsically low luminosity the compan- 
ions must also be faint, and are frequently rare 
or interesting objects. As an example, WDs 
are idea l targets for the direct detection of 
planets dDebes et al]|2005l; IFarihi et alJ boQi ' 



Hogan et al.1 120091 : I.Mullallv et al.1 l2009h and 
brown dwarf stars (e.g. IFarihi et al.ll2005l ). 

The Sloan D igital Sky Survey (SDSS; 
York et al.l lioooh has increased the numb 
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of spectroscopically identified WDs from two 
thousand to a few tens of thousands. This 
has, in turn, allowed follow-up surveys of spe- 
cific types of WD sy s tems, from pulsat ors (e.g. 



Mullally et al.ll2005l ; iNitta et al.1 120091). to ex- 



tremely low mass WDs (iKilic et al.ll2007al l to 
binaries involving main-sequence stars (e.g. 



Silvestri et all l200fil : iHeller et al.l l2009lh and 



binar ies with neutron stars (jAgiieros et al 
20091 k 



Binary WDs hold the solution to an endur- 
ing problem in astrophysics; the progenitors 
of Type la Supernova (SNIa). The origin of 
SNIa is of great interest given their role in 
galactic chemical evolution and determining 
the nature of dark energy. If a WD accretes 
enough material that its mass approaches the 
Chandrasekhar limit (~1.4Mq), the star can 
no longer be supported by electron degeneracy 
pressure and explodes as a supernova. This 
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scenario explains the lack of observed hydro- 
gen in the spectra of SNIa, as well as the strik- 
ing similarities in the lightcurves and spectra. 
However, the source of the accreted material, 
remains a subject of active debate. 

In the do uble degenerate scenari o for SNIa 
proge nitors (jlben k, Tutukovl 11984 ; IWebbink 
[198J), two WDs in a tight binary in-spiral 
due to the emission of gravitational radia- 
tion. If the total mass is near the Chan- 
drasekhar mass, the merger results in a su- 
pernova. While theoretically appealing, it 
is not clear that nature favors this met hod. 
The SPY survey (jNapiwotzki et al.ll200lh . a 
high precision radial velocity survey of over a 
thousand WDs, failed to find any candidates 
with periods short enough to merge within 
the lifetime of the Galaxy, a nd masses large 



enough to explode as SNIa (jNelemans et al 



2005 



Napiwotzki et al.ll2007h . 



The SDSS offers an opportunity to build 
on the SPY sur vey with a much larger sam- 



ple of stars. Kleinman et al. ( 20041 ) and 



Eisenstein et al.l 12006) meticulously collected 



and classified the spectra of nearly 10,000 
WD and sub-dwarf stars, of which approx- 
imately 8,000 were single stars with hydro- 
gen or helium atmospheres (DAs and DBs 
respectively). The spectra were obtained as 
a series of 3 or more 15 minute exposures 



usua lly taken consecutively (jAbazaiian et al 



20091 ). which makes it possible to identify mas- 



sive companions with orbital periods of a few 
hours or less and radial velocity amplitudes 
>170km.s _1 . The low luminosity of the WD 
means that any fainter companion must be 
a degenerate object (WD, brown dwarf, neu- 
tron star, etc.) or a very late M star; any 
other object would be more luminous than 
the white dwarf. 

SWARMS (the Sloan White dwArf Radial 



veloc ity data Mining Survey, iBadenes et al 



20091 ) exploits these individual exposures to 



double degenerate white dwarf (DDWD) sys- 
tems. Our survey is complementary to the 
SPY survey in that it has a lower radial ve- 
locity sensitivity, but is still sensitive to white 
dwarf companions for many thousands of ob- 
jects. In this paper we present two binary sys- 
tems, SDSS J143633.29+501026.8 and SDSS 
J105353.89+520031.0, with periods of 1.1 and 
0.96 hours respectively. These systems consti- 
tute the shortest period non-interacting dou- 
ble degenerate binaries yet found, and are sig- 
nificantly shorter than the 1.46 hour period 
of the previous rec ord holder, WD0957— 666 



( Moran et al. 19971 ). As there are no visi- 



mine the SDSS spectroscopic database for 



ble absorption lines from the companions our 
mass estimates are only lower bounds, but in 
each case the companion is most likely another 
WD. 

2. Observations and Reductions 

We identified SDSS 1436 (ff=18.2, plate- 
mjd-fiber=1046-52460-594) and SDSS 1053 
(5=18.9, 1010-52649-12) as hydrogen atmo- 
sphere (DA) WDs potentially possessing short 
period companions as part of an on-going sur- 
vey for DDWDs. Although we see radial 
velocity variations between different expo- 
sures, no companion is visible in the spec- 
trum. To confirm these systems as binaries, 
and to measure the orbital parameters, we 
observed both stars with the Dual Imaging 
Spectrograph (DIS) with the 3.5m telescope 
at Apache Point Observatory over 4 nights 
between 2009-02-05 and 2009-02-14. We used 
the B1200 grating with a 1.5" slit for a dis- 
persion of 0.62 A per pixel and a resolution of 
1.8 A FWHM. Each exposure was 10 minutes 
in duration and bracketed by an exposure of 
the Helium, Neon and Argon arc lamps. We 
took several exposures of the spectrophoto- 
metric standard Feige 67 each night to flux 
calibrate our spectra. 

We performed an optimal spectroscopic re- 
duction of each spectral image and flux stan- 
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Fig. 1. — De-shifted and co-added spectrum 
of SDSS 1436 based on observations made at 
APO. The solid red line is the best fit model 
used to estimate the temperature and gravity. 
The average S/N per pixel in this spectrum is 
30, and «7 in each individual spectra 



dard using standard long-slit IRAF routines. 
To optimize the wavelength solution we trace 
the arc lamp spectra with the same trace 
used to extract the corresponding WD spec- 
tra. A prominent Mercury emission line at 
4358 A (courteously provided by the residents 
of White Sands, NM) confirmed our wave- 
length offset to better than our resolution. 
The final flux calibrated spectra obtained 
from the blue CCD spans 3790-5020 A. We 
show an average spectrum from all 4 nights 
for SDSS 1436 in Figure E 

The pressure broadened lines of a DA 
white dwarf are well modeled as a combi- 
nation of a Gaussian co re with Lorentz wings 
(jThompson et al.ll2004l ). We convert the cen- 



troid shift of H 7 to a velocity, and fit a sine 
curve with constant offset to the radial veloc- 
ity time-series. Radial velocity curves for Ha 
and H,5 give similar results, but the accuracy 
obtained by fitting H 7 alone is sufficient for 



our purposes. We show the best fit folded ra- 
dial velocity curves in Figures [2] & [3] and the 
best fit parameters in Table [TJ Given the more 
sparse sampling of SDSS 1053, our period es- 
timate is less certain than for SDSS 1436, but 
our observations span nearly 2 orbits and our 
uncertainty estimate is only 36s. 

The residuals of the fit to SDSS 1436 show a 
linear trend with phase. This trend is seen for 
fits to and as well. Examination of the 
unfolded lightcurve confirms this trend is in- 
deed a function of phase, not of time, and can 
not be explained by some drift in our instru- 
mental calibration. Similarly, a third body in 
the system on a longer period orbit would only 
produce a trend in the unfolded data. Fitting 
an eccentric orbit reduces the peak to peak 
amplitude of the residual trend by 100 km.s -1 
but does not eliminate it. Because the eccen- 
tric orbit fit is not significantly better, and 
because we have difficulty imagining a sce- 
nario in which a system that has undergone 
two common envelope evolution phases could 
emerge with an eccentric orbit, we show only 
the circular fit in Figure [2j 

2.1. Temperature and Gravity 



Th e published temperatu re of each star 



from 
fit to 



Eisenstein et al 



(|2006h comes from a 
the average of three separate expo- 
sures spanning a total of 50 minutes. As 
this is a significant fraction of the orbital pe- 
riod we were concerned that the fit may have 
been biased by combining spectra with dif- 
ferent radial velocities. Using our best fit 
radial velocity curve, we deshifted and co- 
added each of our 10 minute exposures to 
produce a high signal-to-noise spectrum. We 
then fit this spe ctrum to the sam e grid of DA 
models u sed by Eisenstein et al. ( 20061 ) (up- 
dated by Koester et al. 20091 . and kindly pro- 
vided by the author). We linearly interpo- 
lated the model spectra to produce a finer 
grid of AT efr =10K and A log #=0.02. We 
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Phase (n) 

Fig. 2. — Radial velocity curve for SDSS 1436 
folded at the best fit period. The solid line is 
the best fit sine curve to the data. The resid- 
uals of the fit are shown in the lower panel. 
The apparent linear trend in the residuals is 
discussed in the text. 



fit each model to the entire spectrum from 
3800-5000 A using a least squares minimiza- 
tion algorithm, allowing the fit to vary by a 
high order polynomial in a similar manner to 
Eisenstein et al. I tod ). We find best fit pa- 
rameters of 17120 K, log #=6.60 for SDSS 1436 
and 16150 K, 6.35 fo r SDSS 1053, consisten t 
with the estimate of Eisenstein et al. (2006) 
who finds (T e ff, logg) of (16933, 6.58) and 
(15399, 6.28) respectively. 

For each object we combine spectra taken 
close to the minima and maxima of the ve- 
locity curve. We find no evidence of spec- 
tral features in these combined spectra and 
are confident that flux fro m the compan 



i on is not biasing our fit. iFontaine et al 
(|2003h noted that temperature and gravity 
estimates of WDs from independent spectra 
using identical reductions and identical at- 
mosphere models often disagree significantly 
more than the quoted uncertainties. Follow- 
ing their approach, we adopt uncertainties of 




Fig. 3.- 
SDSS 1053. 



Phase (n) 

Same as Figure El but for 



200 K and 0.05 for our fits, which are more 
conservative than the values returned by the 
fitting method. We caution that these uncer- 
tainties are internal to our fitting, and do not 
attempt to a ddress limitation s of th e models. 
For example, Tremblav et al. ( 20091 ) recently 
introduced an improved treatment of Stark 
broadening which systematically increases the 
best fit gravity by 0.2 dex in this temperature 
and gravity range. 



Kilic et al.l (|2007al ) independently observed 



and fit the spectra of both stars (as part of 
a search for companions to low mass WDs) 
and obtained similar results for the gravity 
(log 5= 6.59 and 6.40), but higher temper- 
atures (T e g= 18339 and 18325). Given the 
close agreement in measured gravity between 
the three measurements, the small discrep- 
ancy in temperatures do not materially effect 
the stellar masses we estimate in Section [3] 

We simulated the effect of changing radial 
velocities over the course of a 10 minute ex- 
posure to estimate the effect on the best fit 
temperature and gravity. Using our best fit 
radial velocity curve for SDSS 1436, we coad- 



ded a series of appropriately Doppler shifted 
model spectra, and fit the result in a man- 
ner similar to our data. The largest discrep- 
ancy occurs when the star is traveling per- 
pendicular to the line of sight, where the core 
of the line appears smoothed. This blurred 
spectrum is preferentially fit by 500 K hotter 
model with a shallower line core, but the best 
fit gravity, which is most important for mea- 
suring the mass, remains unchanged. 

3. Discussion 

Comparing our best fit temperat ure and 



gravit y to the WD evolution models of lSerenelli et al 



(j2002h we estimate masses of 0.23(01) and 



0.21(01) Mq. These models are created by re- 
moving mass from a 1 Mq model at appropri- 
ate times during red giant branch evolution, 
and incorporate chemical diffusion, a nearly 
pure He core (with metallicity, Z=0.001) and 
thick H layer. The estimated masses are close 
to the minimum known white dwarf mass of 
0.17 M, 



200Sh . 



> (Kilic et al. 2007a: 


Kawka & Vennes 


Moroni & Straniero ( 


2009) estimates 



that a WD must have a mass of at least 
0.33 Mq to have a carbon-oxygen core. Nei- 
ther object approaches this mass, and are 
composed almost entirely of helium with a 
thin hydrogen atmosphere. 

Accord ing to the initial-final mass rela 



tion (e.g. Kalirai et al. 20081 ; Williams et al 



20091 ). only isolated WDs with masses greater 
than 0.47 Mq have had time to evolve off 
the main-sequence within the lifetime of the 
Universe. Although it has been argued that 
high metalli city progenitors ca n produce lower 



mass WDs (jKilic et al.ll2007d ). that scenario 



is unnecessary for these systems. Instead, the 
fact that these two systems are known to be 
binaries suggests that growth of these lower 
mass WDs was truncated by a common enve- 
lope pha se and evolve d through the sub-dwarf 
channel (|Heberll2009l ). 
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Fig. 4. — Distribution of known DDWD sys- 
tems. The filled circles indicate the mini- 
mum mass of the two new systems discussed 
in this paper. The notch on the arrow gives 
the total mass assuming an inclination an- 
gle, i, of 60°, and the tip of the arrow 
shows i = 40°. Square symbols indicate 
previously known double-lined (DL) binaries 
(where the total mass is known), while tri- 
angles indicate single-lined (SL) systems and 
are only a lower bound on the total system 
mass. SPSS J0917 09.55+46382 1 .8 (op en tri- 



angle) taken from Kilic et al. ( 2007bh. an d 



all other systems from lNelemans et alJ Q. 
The horizontal dashed line indicates the 
Chandrasekhar mass, while the curved line 
shows the period for which the merger time 
is equal to the age of the Universe. 



3.1. Nature of the Companions 

Solving for the Keplerian equations of 
motion for SDSS 1436 gives a mass for the 
companion of (0.57(04)/sin l) Mq, where t 
is the inclination of the orbit to the line of 
sight, consistent with a carbon-oxygen core 
WD. The companion to SDSS 1053 is at least 
0.31(02) Mq. Although these are minimum 
masses, and are consistent with a wide range 
of astrophysical objects, we argue that the 
companions are most likely also WDs. 
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Table 1 
Measured Parameters 





SDSS 1436 


SDSS 1053 


Period (hrs) 


1.15238(14) 


0.960(10) 


Amplitude (km.s -1 ) 


388(21) 


310(14) 


Seperation (Rq) 


0.4789(75) 


0.2924(51) 


Teff(K) 


17120(200) 


16150(200) 


log g 


6.60(05) 


6.35(05) 


Massi (M ) 


0.23(01) 


0.22(01) 


Mass 2 /sini (M ) 


0.57(04) 


0.31(02) 


Merge time (Myr) 


<102 


<104 



Main sequence stars can be ruled out on lu- 
minosity grounds. For SDSS 1436 (SDSS 1053), 
the minimum mass of the companion cor- 
responds to a spectral type of K8 (Ml) 
dHabets k Heintzjll98lh . which has an abso- 
lute % magnitude of 7.2 (8.5) dBilir et alJ feoOQ: 



Hawlev et al.l |20Q2| ) . considerably brighter 
than the observed WD (i = 9.1 (9.4), 
Holberg Bergeron! [20061 k The SDSS spec- 
trum of either object shows no evidence of any 
cool companion at red wavelengths, ruling out 
the possibility of a main-sequence companion. 
A similar argument applies to red giant stars 
and other, higher luminosity objects. 

If the inclination angle, t < 24° (13°), the 
companion mass is greater than the Chan- 
drasekhar mass and the companion must be 
a neutron star (NS) or a black hole. Ap- 
proximately 45 WD-NS binaries are known, 
and the mass distribution of WDs in such 
binaries is much wider than for isolated sys- 
tems, admitting both high and low mass WDs 



van Kerkwiik etaJ] 120051 k lAgiieros et al 



2009J.) looked at both SDSS 1436 and SDSS 1053 



during an 820 Mhz radio survey for pulsar 
companions to WDs but did not detect any 
signal. These observations do not exclude the 
possibility of a pulsar companion, not only 



because the orientation of the pulsar beam 
may not be along the line of sight, but also 
because their analysis restricted their sensi- 
tivity to orbital periods greater than 8 hours. 

Interacting WD-NS binaries are known as 
ultra compact X-ray bin aries (UCXBs; see 
Nelemans Sz Jonker 20061 . for a review). In 



these systems, the orbital separation is so 
small that a WD overfills its Roche lobe 
and donates material onto the surface of 
the neutron star via an accretion disk, emit- 
ting X-rays in the process. The longest pe- 
riod UXCBs have periods of 50-55 minutes 
( Nelemans fe Jonker 20061 ). entirely consis- 
tent with the periods of the systems under 
scrutiny. If the companions were NSs, they 
would almost certainly be interacting. How- 
ever, the presence of hydrogen in the atmo- 
sphere of the visible WD in both systems 
means that the systems are not interacting. 
In double degenerate systems, mass trans- 
fers from the lower to the higher mass star. 
In both systems, the higher mass star is the 
invisible companion. If the visible star was 
losing mass, the thin hydrogen layer would be 
quickly stripped, exposing the underlying he- 
lium core. Because both stars still have their 
hydrogen layers we can conclude that mass 
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transfer has not yet started. 

3.2. Consequences of a Merger 

Non-interacting DDWD companions are 
therefore the only possible objects consistent 
with the available evidence. Non-interacting 
systems in short period orbits lose orbital en- 
ergy in the form of gravitational radiation and 
will eventually merge. Using the e quation for 
angul ar momentum loss given by iPaczvhski 

ESS, 

we estimate these systems will merge 
in less than 102 and 104 Myr respectively. If 
the orbits are inclined by 45° to the line of 
sight, the merger time decreases to < 79 Myr 



Guerrero et al.l (|2004h used smoothed par- 



ticle hydrodynamic simulations to predict the 
consequences of the merger of WDs of differ- 
ent masses. For the merger of two 0.4 M Q he- 
lium core WDs (analogous to the minimum 
mass configuration of the SDSS 1053 system), 
they found no thermonuclear flash, and no 
mass loss. The case of a 0.4 M WD merging 
with a 0.6 M carbon-oxygen core star (simi- 
lar to SDSS 1436), some carbon is burned into 
oxygen, but there is no thermonuclear run- 
away and no supernova. 

Unless the unseen companions have masses 
> 1.2 Mq it seems unlikely that either sys- 
tem will produce a SNIa. WDs with masses 
greater than 1.2 Mq do exist, with the most 
massive WD found to date being 1.33 M 



(jKepler et al.1 120071 . assuming an oxygen-neon 
core). However, these stars are rare, and are 
consistent with only a small range of inclina- 
tions angles for these systems. 

The merger of a carbon-oxygen core WD 
with a helium WD most likely produces an 
extreme helium star (ISaio &: Jeffervl 120020 or 
an R CrB star (jWebbinklll984l : Iciayton et al 



20071 ). The merger of two helium core WDs 
(an option only for SDSS 1053) is one evo- 
lution ary pathway to produce sdO subdwarf 
stars (jHeberl l200a ). Regardless, the merger 



remnant will eventually cool t o become a sin- 
gle WD. iLiebert et ail ifcoOfih compared the 
space density of high mass WDs with that ex- 
pected from a Salpeter initial mass function 
and a single burst stellar population and con- 
cluded that 80% of high mass WDs were cre- 
ated by the merger of lower mass stars. The 
merger of a carbon-oxygen WD with a He core 



one ha s been suggested by lGarcia-Berro et al 



(120071 ) as the origin of hot debris disks around 
massive WDs. 

However, analysis of the extremely low lu- 
minosity, and calcium rich, t ype l b supernova 



SN2008E by IPerets et al.1 (120091 ) concluded 



that only 0.3 M of material was ejected, and 
the pattern of elemental abundances (high cal- 
cium abundances, but low sulpher) was best 
explained by helium fraction >0.5 in the ini- 
tial composition. It is conceivable that the 
progenitor of this explosion involved the dis- 
ruption of a helium core WD in a DDWD sys- 
tem. 

4. Conclusion 

We report on the detection of the two clos- 
est non-contact white dwarf binaries known. 
These systems were detected by mining the 
spectroscopic database of the SDSS, and 
followed-up with time resolved optical spec- 
troscopy. We argue that the companions must 
also be WDs: Main-sequence stars of the req- 
uisite mass are more luminous than the pri- 
mary WDs, and a neutron star or black hole in 
such close proximity would have stripped off 
the thin outer layers of the primaries. With 
periods of about an hour, these systems will 
merge in less than 100 Myr and probably pro- 
duce a high mass WD. 
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